JOURNAL DE PHYSIQUE IV 

Colloque C3, supplement au Journal de Physique H, Volume 3, aout 1993 




Tractable chemical models for CYD of silicon and carbon 

E. BLANQUET^ 1 ) and SA. GOKOGLU 

NASA Lewis Research Center, 21000 Brookpark Rd, Cleveland OH 44135, U.SA. 


ABSTRACT 

Tractable chemical models are validated for the CVD of silicon and carbon . Dilute silane 
(SiH4) and methane (CH4) in hydrogen are chosen as gaseous precursors. The chemical 
mechanism for each systems Si and C is deliberately reduced to three reactions in the models: 
one in the gas phase and two at the surface. The axial-flow CVD reactor utilized in this study 
has well-characterized flow and thermal fields and provides variable deposition rates in the axial 
direction. Comparisons between the experimental and calculated deposition rates are made at 
different pressures and temperatures. 

INTRODUCTION 

Chemical Vapor Deposition (CVD) is a widely used process in modem technology for the production of 
many electronic, optical and structural materials. Despite its increasing importance in the industry, the 
scientific understanding of the phenomena involved in CVD is still rather insufficient. A complete 
analysis of this process should include gas phase and surface chemical reacdon kinetics, multicomponent 
mass transport, fluids mechanics and thermodynamics [1]. Many different numerical models have been 
developed with broad applications to the CVD process. The availability of elaborate computational fluid 
dynamics (CFD) codes allows accurate predictions of velocity, thermal and concentration fields in three- 
dimensional or complex geometry reactors. It appears that the correct chemical description of the global 
approach is the weakest link of the chain today. This arises mainly from the difficulty in either obtaining 
reliable gas phase and surface kinetic information and/or proposing the accurate chemical mechanisms 
evolving in the CVD reactor. Chemical models treating hundreds of .gas phase and surface reactions 
affect the degree of accessibility and reliability of such models [2], 

APPROACH 

■fhis paper presents a study involved in a systematic experimental and numerical approach -described 
elsewhere [3] - in order to propose tractable chemical models for the system Si -C -H. In our systematic 
approach, we first study the CVD of the individual elemental constituents of SiC, namely, silicon and 
carbon. Dilute silane (SiH^) and methane (CH4) in hydrogen are chosen as gaseous precursors for CVD 
•of silicon and carbon in an effort to provide analogous chemical mechanisms and reaction pathways. The 
chemical mechanisms for the two systems are deliberately reduced to three reactions: one gaseous and 
two surface reactions as schematically described on figure 1. An imaginary intermediate species A m Hn is 
included in the mechanism path. The identity and subsequent reactions of the imaginary intermediate 
species A m H„ will certainly affect the magnitude of the effective rate constants k* and T|'. These 
effective species and rate constants created by our "lumping" strategy are expected to take into account 
the overall actual reactions. 


(1) National Research Council Resident Research Associate. Present address: LTPCM-ENSEEG, 
Domaine Universitaire, BP. 75, 38402 Saint-Martin d’H'eres cedes, France 
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AH4 


k* : effective gas phase reaction rate 
constant consuming AH 4 

Intermediate species A mHn 


TV 


T)' : effective surface reaction rate constant 
for the conversion of A mHn to solid A 


ti : effective surface reaction rate constant 
for the conversion of AH 4 to solid A 


surface species A = Si or C 


Figure 1. Chemical mechanism for CVD of AH 4 . 

The finite-volume based computational code FLUENT 3.03 [4] is adopted to simulate the CVD reactor in 
a two-dimensional flow geometry. The reactor has well-characterized flow and thermal fields and is 
designed for variable deposition rates. It is a hot-wall, axial flow reactor schematically shown in figure 
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Figure 2. Schematic of the hot-wall, axial-flow horizontal CVD reactor. 


It provides simultaneous variation in gas phase concentration and temperature. A thermocouple alumina 
tube (OD= 6.3 10 " 3 m) is installed in the center of the horizontal heated alumina reactor 
(ID= 2.54 10-2 m) and is used as substrate holder. Temperature profiles on the surface of the external 
alumina tube and inside the substrate holder, total pressure and inlet gas mass flows are measured. 
Deposition is done on Saphikon™ and Textron™ fibers, cemented on top of the centered alumina tube, 
as shown in figure 2. Prior to deposition, fibers are ultrasonically cleaned and then blown dry. Deposit 
film thickness is measured along the fiber by an optical method (Laser Mike™) using the shadow of a 
laser beam and verified by Scanning Electron Microscopy. For silicon CVD, two temperatures ("1000 
and 1200°C") corresponding to the maximum of temperature profile of the furnace are selected. The 
concentration of silane is 0.3% by mole in hydrogen. The temperature "1400°C" and methane ! 
concentration of 5% by mole in hydrogen are chosen for carbon CVD study. In both cases, total 
pressure is either 9.3D10 4 Pa (700 Torr) -labeled AP- or 6650 Pa (50 Torr) -labeled LP. Total flow i 

rate is kept constant at 1 slm in all runs. For each conditions, at least two runs have been performed and 
uncertainties related to thickness measurement have been considered. 







Temperature (*C) 
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RESULTS AND DISCUSSION 

For the numerical calculations, we use the experimentally measured temperature profiles on the external 
wall and inside the substrate holder for each total pressure and gas flow rate. Geometry of the reactor is 
treated as axisymmetric. Because methane and silane are dilute in hydrogen, thermophysical and 
transport properties are calculated by treating the mixture as pure hydrogen. Reactions are considered on 
both external tube wall and substrate holder surfaces. 

Temperature profiles at the surface of the substrate holder calculated from the experimental values for the 
"1000 and 1400°C" cases are presented with the corresponding experimental deposition rates in figures 
3a and 3b, respectively. The temperature profiles are essentially independent of total pressure. 


Silicon LP "1000“C" Carbon AP "140CTC" 



Figure 3. Calculated temperature profile and experimental deposition rate (not in scale) along the 
substrate, a) Silicon "1000°C", b) Carbon "1400°C". 

For the case of carbon, there is no deposition when temperature is below 1300°C (see arrows in figure 
3). For the case of silicon, deposition rate increases at around 750°C, and drops because of silane 
depletion. 

For carbon deposition, the considered gas phase and surface reactions are given as follows: 

Gas phase: CH4 --> An,H n (Rl) Surface: CH4 — > C (R2) and A m H n — > C (R3). 

Reaction 1 (Rl) is an irreversible reaction. Different values for the reaction Rl rate constant, listed in 
table 1, selected from the NIST Chemical Kinetics Database [5] and the literature review proposed by 
( Warnatz [8] have been considered. Surface kinetics of reaction R2 and R3 are estimated using different 
sticking coefficients 7 in the equation: T)' = ( 7 ) ( RT / 2it m) 1 ^ 2 - where m is the species molecular 
weight. Methane sticking coefficient has been reported to be equal to 5 E-5 [9] and lower [10]. 


Reference 

> 

3 

M 

A ( s' 1 ) 

P 

E (J/kgmol) 

Klotz [6], NIST 

Unidentified 

n 2 

IEI3 

0 

4.28E8 

Klotz [7], NIST 

Unidentified 

h 2 

1E9 

0 

2.43E8 

Wamatz [8] 

5H3 

forward reaction 

- 

1E15 

0 

4.20E8 


Table 1 . Reported rate constant k* for methane dissociation CH4 --> A m Hn , in diluent M 
k*=A TP exp(-E/RT). 
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During the process of parameter optimization, methane sticking coefficient has been deliberately reduced 
to IE-10, which is consistent with Makarov and Pechik results [10]. The predicted values of C 
deposition using the values proposed by Klotz [7] pertaining to a hydrogen bath are much larger than the 
experimental ones, therefore the values proposed by Klotz [6] and Wamatz [8] have been used. Figures 
4a and 4b show the optimized predicted and experimental carbon deposition rates along the substrate 
under the conditions "1400 P C” - LP and "1400°C 1 ' - AP, respectively. 


Carbon LP " 1400°C" Carbon AP " 1400°C" 



Figure 4. Predicted and experimental carbon deposition rates along the substrate under the conditions: 
a) "1400°C", LP b) 'T400 0 C'', AP. The error bars are representative of experimental uncertainties. 

By varying A (Reaction Rl) and y (Reaction R3) the shape of the calculated deposition peak changes: A 
and y have an effect on the maximum of the peak and its position, respectively. In both cases, AP and 
LP, the optimized A value is 2 orders of magnitude lower than the value recommended by Wamatz [8], 
which might be explained with two considerations: a) the value recommended by Wamatz is the forward 
reaction rate whereas our "lumped" value corresponds effectively to both forward and reverse reaction 
rates and b) the working pressure (50 or 700 Torr) belongs to the domain of reaction rate "fall-off' 
(decrease of the rate constant with pressure) for both forward and reverse reactions, as shown in figure 
5. 
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Figure 5. Fall-off curves for a) CH3 + H — > CH4 and b) CH4 «> CH3 + H. Taken from [8] and [11], 
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I For silicon deposition, the considered gas phase and surface reactions are given as follows: Gas phase: 
■ S1H4 < — > S1H2 + H2 (Rl) Surface: SiH} -> Si +2 H2 (R2) and Sffife — > Si + H2 (R3). Reaction 1 
i (Rl) is a reversible reaction. Two different values for the reaction Rl kinetics have been selected: one - 
j labeled, "Model A"- is taken from Moffat and Jensen [12], (kf = 6.1 E28 T*5 exp(-2.461 E8/RT) s _1 and 
k r = 5.28 E21 T -4 - 44 exp(-1.427 E7/RT) m 3 /kgmols , preArrhenius factor in J/kgmol), the second 
labeled "Model B" come from calculations made for higher temperature silane dissociation [13], (kf = 
4.07 E29 T- 4 - 966 exp(-2.505 E8/RT) s' 1 and k r = 6.65 E22 T -4 - 041 exp(-1.3794 E7/RT) m 3 /kgmols 
preArrhenius factor in J/kgmol). For both models A and B, the surface kinetics of reaction R2 is 
calculated from the sticking coefficient of silane given by Coltrin et al. [14] (y= 5.37 E-2 exp(-9400/T)). 
For reaction R3 the sticking coefficient for SiH2 is taken as unity. Figures 6a and 6b show the predicted 
and experimental Si deposition rates under the conditions "1000°C" LP (6a) and AP (6b), respectively. 


Silicon LP "lOOO’C" Silicon AP "1000”C" 



Figure 6. Predicted and experimental silicon deposition rates along the substrate under the conditions: 
a) "1000°C", LP, b) "1000°C", AP. The error bars are representative of experimental uncertainties. 

In both cases Model A provides a good agreement with experimental results for the shape, position and 
magnitude of the deposition profile. Calculations performed with decreasing the temperature by 25° 
(which is an acceptable error for temperature measurement) indicate a better agreement as shown in figure 


Silicon LP "1Q00°C" 



Figure 7: Predicted and experimental silicon deposition rates along the substrate under the conditions: 
"1000 o C", LP. Model A T^TOOCFC" and Model A T=''975°C". 
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• However, for the higher temperature "1200 °C" runs -displayed in figures 8a (LP) and 8b (AP) it 
appears that Model B -as expected- is predicting the magnitude of the deposition peak better than Model. 
A but fails for its position. 


Silicon LP "1200°C" Silicon AP "120CTC" 



Distance From Inlet (m) Distance From Inlet (m) 


Figure 8. Predicted and experimental silicon deposition rates along the substrate under the conditions: 
a) "1200°C", LP, b) "1200°C", AP 


Similar to the "1000°C" cases, the effect of some experimental parameters has been investigated: 
a) temperature decrease of 200°C (figure 9a), b) silane concentration reduced from 0.30% to 0.17% 
(figure 9b) and c) misplaced -by three cm- temperature profile measurement (figure 10). It is shown that 
these acceptable changes in boundary conditions would promote either of the two models. 


Silicon LP "1200°C” Silicon LP "1200“C" 



Distance From Inlet (m) Distance From Inlet (m) 


Figure 9. Predicted and experimental silicon deposition rates along the substrate under the conditions: 
”1200°C", LP. a) Model A change in T, b) Model A change in silane concentration. 
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Figure 10. Predicted and experimental silicon deposition rates along the substrate under the conditions: 
"1200°C", LP. Model B temperature profile misplaced by 3 cm. 

j *» •. 

SUMMARY AND CONCLUSION 

j 

j A systematic experimental and numerical approach is designed to develop tractable chemical models for 
CVD applications. The model is assumed to consist of lumped chemical reactions, one in gas phase and 
two at the surface. Initial testing of the approach is carried out by using the chemical kinetic parameters 
currently available in the literature for silicon and carbon deposition from silane and methane dilute in 
hydrogen. Comparison of the model performances with our experiments show that adjustments are 
■ necessary to the kinetic parameters for better predictions. Optimized tractable models consistent with 
literature data are proposed. 
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